PACS numbers: 98.80. Jk, 96.12.Fc 



Singularity problem in f(R) model with non-minimal coupling 

Qiang Xu^Q and Sheng-yu Tan 1 ^ 

1 Department of Physics, and State Key Laboratory of Nuclear Physics and Technology, 
Peking University, Beijing 100871, P.R. China 

We consider the non-minimal coupling between matter and the geometry in the f(R) theory. In the 
new theory which we established, a new scalar tjj has been defined and we give it a certain stability 
condition. We intend to take a closer look at the dark energy oscillating behavior in the de-Sitter 
universe and the matter era, from which we derive the oscillating frequency, and the oscillating 
condition. More importantly, we present the condition of coupling form that the singularity can be 
solved. We discuss several specific coupling forms, and find logarithmic coupling with an oscillating 

■ period AT ~ Az in the matter era z > 4, can improve singularity in the early universe. The result of 
y—\ ' numerical calculation verifies our theoretic calculation about the oscillating frequency. Considering 

two toy models, we find the cosmic evolution in the coupling model is nearly the same as that in 
£S) ■ the normal f(R) theory when Ina > 4. We also discuss the local tests of the non-minimal coupling 

" f(R) model, and show the constraint on the coupling form. 

c/r 
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"^J 1 - I. INTRODUCTION 

cr 

The rapid development of observational cosmology starting from 1990s shows that the expansion of our universe 
bJ[) in the present epoch is accelerating. The most popular theory to explain the accelerating universe is dark energy 
(DE). One possible dark energy candidate accepted by most of researchers is the cosmological constant. However, the 
^{■^ , cosmological constant that remains a mystery, can not be explained clearly in any known theory. If the cosmological 
constant originates from the vacuum energy in quantum field theory, as many people believe, its energy scale is too 
Q\ \ large to be compatible with the observed dark energy density Moreover, the observation indicates that the dark 
t-H ■ energy equation of the state may cross the phantom divide line lj = — 1. This suggests that the cosmological constant 
may not be the only candidate for dark energy. There also exist some other dark energy models, ranging 
from quintessence, phantom, quintom to chaplygin gas models. For the nice reviews on dark energy, see An 
0^ ' alternative scenario for dark energy is infrared(IR) modified gravity. Among many IR modified gravity models, f(R) 
gravity is of particular interest. One important feature in f{R) gravity is the intrinsic existence of an extra dynamical 
scalar degree of freedom, besides the massless graviton. Therefore it is possible to study both the early-time inflation 
and late-time acceleration of the universe in the framework of f(R) gravity, without introducing a scalar field by hand. 

■ More interestingly, the effective equation of state could be smaller than —1 in f(R) dark energy models, indicating 
the scalar behaves like a phantom in the Jordan frame. The model with a Lagrangian density f(R) = R — a/R n 
(a > 0, n > 0) was proposed for dark energy 

d-d- 

However this model is plagued by matter instability [l2l . [iH and 
difficulty to match local gravity constraints. Later on, researchers have proposed many viable models, seeing [l 4M22j j . 
For nice reviews on f(R) theories, see [23l |24|. 

Nonetheless, it has been pointed out in Refs. [U I25rl3f3 | that viable f(R) models generally suffer from a singularity 
problem. For most of the f(R) models, in order to evade the local gravity tests, Rf" < 1 is a common feature. 
This leads that the oscillating frequency F(z) ~ 1/^/Rf" becomes very large. In this case, the Ricci scalar, Hubble 
parameter and EOS parameter of DE have intense oscillations in the high redshift region. One possible way proposed 
by [3l[ to solve this problem is adding a R 2 /M 2 term to the f(R) function. However, according to the constraint 
of f(R) inflation, M ~ 10~ 6 M p i. If we consider R ~ A, then the adding term contributes A/M^ — 10~ 121 to the 
oscillating frequency. So, this term may has some contribution in high red-shift region such as z > 10 40 , but in the 
matter era z > 4, such term is of little help. Even recently, this problem is till be discussed in [32|. 

A generalization of the f(R) theories was proposed in [33l-l35j firstly by including the theory an explicit function of 
the Ricci scalar R with the matter Lagrangian density L m . As a result of the coupling, shown in (36|, the equation of 
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motion which is the non-geodesic, and an extra force orthogonal to the four-velocity, arise. The implications of the 
non-minimal coupling on the stellar equilibrium were discussed in 37], where the constraints on the coupling were 
obtained. The equivalence between a scalar theory and the model with the non-minimal coupling was considered in 
[38l l39l ] , where the authors showed that the non-minimal coupling f(R) theory corresponds to a two-field scalar theory. 
Especially in the non-minimal coupling model, the matter part was extended to a arbitrary function of the Lagrangian 
density of the matter in [40l |. Later, the coupling model was analyzed in [4l| to study the matter perturbation and 
gave a possible way to prevent the f(R) theory conflicting with the galaxies matter spectrum tests, in [52] to study the 
accelerated expansion of the universe, and in [43j to discuss the reheating, which gives a constraint A/i?i < 10 -104 
about the coupling form f 2 = 1 + (R/Ri) n . 

In this work, we study the oscillating behavior of the dark energy in the non-minimal coupling f(R) theory. We 
consider the non-pressure dust as the main contribution of the matter Lagrangian and the coupling form is arbitrary 
function of Ricci scalar. So the energy-momentum tensor of the matter is generally conserved and the matter density 
is proportional to 1/a 3 . Similar to the f(R) theories, the non-minimal coupling f(R) model also has a scalar freedom 
-0 which consists of the f(R) part /{, and the coupling part, pf 2 . Through adding the coupling part, the oscillating 
frequency F(z) ~ 1/ \/Rf" can be modified, so that the singularity problem may be solved. 

This paper is organized as follows. In Section II, we present the definitions in the non-minimal coupling f(R) model 
and derive the equation of motion. Then we study the stability condition. In Section III and IV, we solve the EOM 
of dark energy in de-Sitter universe and the matter era theoretically. And we mainly focus on the study of oscillating 
behavior of the EOS parameter of dark energy. In Section V, we show that a logarithmic coupling can solve the 
singularity problem. In Section VI, we give the numerical result, and discuss some toy models. In Section VI, we 
discuss the local tests of the non-minimal coupling f(R) model. Finally, in section VIII, we make the conclusions. 



II. DEFINITIONS AND THE STABILITY CONDITION 



A. definitions and the conservation equations 



The action [33J we consider is 



d 4 x, 



1 

2k 



2 h{R) + h{R)L 



(i) 



The parameter k 2 = 87rG, and we set it to be unit in the following sections. L denotes matter Lagrangian. By varying 
the action with respect to the metric g^, following (3(| SH SHwith some modifications, we get the modified Einstein 
equation 



{F x + 2LF 2 )R liV - l -h9^ + (g^a - VpVOCFl + 2LF 2 ) = foT^, 

where F[{R) — f-(R). The matter energy- momentum tensor is defined as 

2 6(^7jL) 



T — — - 



-g 5{g^) 

Using the Bianchi identity, V^G^ — 0, and the identity 

(□V, - V v n)Fi = R^VFt, 
and following (3(| SS SH SE| ; we deduce the following covariant conservation equation 

F 2 



(g^L-T^VR, 



(2) 



(3) 



(4) 



(5) 



which indicates the non-minimal coupling between curvature and matter yields a exchange between matter and 
the geometry. In the absence of the coupling, f 2 (R) = 1, one recovers the covariant conservation of the energy- 
momentum tensor. A simple choice of the matter Lagrangian is L — —p, if we consider the non-pressure dust as the 
main contribution. Considering the FRW metric, T 00 = p and g 00 — — 1, we find the matter density is conserved 



VqT 1 



no 



F 



2 (g 00 L - T 



()0\ 



0, 



(6) 
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which yields the conserved equation of the matter density 

p + 3Hp = 0. (7) 

So, the evolution of the matter density is the same as that in ACDM, p = poa~ 3 . 

In a fiat FRW metric with a scale factor a, according to Eq. @, we get the modified Friedmann equations, 

WH 2 = / 2P+ i(^_/ 1 )_3if^, (8) 

-2^ = / a p+(*-ir|)* (9) 

where tp = Fx + 2LF 2 , H = -, R = 12H 2 + 6H and the dot denotes a derivative with respect to the cosmic time t. 
By rewriting Eq. ([8]), we can define the effective energy density, 

Peff =3H 2 = ±- if 2 p + l(#R - fx) - 3ff^| , (10) 

and the effective pressure, 

P eff = - { 2H + 3H 2 ) = I L\ m - f x ) + ( * + 2H±\ A . (11) 



We can define the dark energy density pd e as pde — Pef / — /?, and we use the definition in the previous papers |14l . |21| , 

y( ^g! = g!_ ( i + ^3. (12) 



Here m 2 is the mass scale m 2 = K 2 p /3 and po is the matter density at the present time. 
The EOS-parameter for dark energy can be expressed by y(z) as 



= -1 + -(1 + z) . (13) 

6 y(z) dz 



When y{z) tends to be a constant, the EOS-parameter tends to be -1. 



B. field equations and the stability condition 

In this section, we discuss the stability in the local gravity. Taking the trace of Eq. ^ , we get 

3n^ + ^-2A(i?)-/ 2 T = 0. (14) 

First, we decompose the quantities R, ip(R) and into the background part with a constant curvature and the 
perturbed part: R = Rq + 5R, ip = ipo + Sip and T = Tq + ST. We consider R close to the mean- field value Rq, and 
the metric still very close to the minkowski case. The linear expansion of Eq. (| 14[) in a time-independent background 
gives 

V 2 SiP - M|<5V = ^-ST, (15) 

where is the mass of the scalar 

^■KtT- 1 )' (16> 

and ip' = F[ + ITqF'^. The stability condition is given by 

0<_^_<i. (17) 
Fx — F2T0 

Note that, to respect the solar system constraints, we must let / 2 — > 1, Fx — > 1 and F 2 —> 0. So this condition (|T7|) 
can be rewritten approximately as 

< Rip' < 1. (18) 

Especially in the late time de-Sitter universe, p = —To — > 0, - > -Pi and V'' ~~ ► then the stability condition p7|) 
recovers to that in the normal f(R) gravity. 
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III. OSCILLATIONS IN THE DE-SITTER UNIVERSE 

The trace equation (fl"4]) of the field can be recast in the form 

□* = ^. (19, 

Similar to the dynamic of normal f(R) gravity, there is a new scalar freedom tp which decides the cosmic evolution, 
and the effective potential has a minimum at 

^ = ~ [i# ~ 2A - / 2 T] = 0. (20) 
dip 3 

In the de-Sitter universe, neglecting the contribution of the matter, the effective potential has a minimum at 

ity - 2A = 0, (21) 
where Rds is a constant. Using Eq. (fT2|) . the Ricci scalar can be expressed as 



R = 2>wf [ 4y(z) - (1 + z)^^- + (1 + zf ) , (22) 



By combining Eq. flSJl with Eq. (|22|). one gets 

d 2 y{z) a dy(z) b , . . . d 



dz 2 1 + z dz (l + z) 2tfw v y (1 + 2) 



= e(l + z) + ■ , ? , (23) 



where 



2 — -0 

o 1 - V» 6F 2 



(25) 
(26) 



d = R ~l 1 (27) 
6H 2 3m 2 ip' ' 

and we have used d/di = —(1 + z)H(z)d/dz. 

When z — > — 1 + , we consider the perturbations around the de-Sitter solution of the dark energy density y(z) — » yo, 
hence 

R^R ds ~ I2m 2 y ~ 12tf 2 S , (28) 

p = 3m 2 (l + z) 3 ^0, (29) 

ip = F t - 2 P F 2 -> Fi, V' -> F{. (30) 

In this case, we can rewrite the coefficients as 

= _3-fc#), (31) 



4(2 -V>) 



(32) 



d = 2 S~y. (34) 



The solution of Eq.® is 



= y + 5(1 + z )*U-°±V(-D 2 -*>) + -1— (1 + z) 3 , (35) 

3a + b + 6 
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where B is a constant depending on the initial condition and yo = d/b. Depending on the sign of the discriminant 
in the square root of Eq. (I3"5j) , there are two possible behaviors for this model. If (a — l) 2 — 4b > 0, the solution 
approaches the de-Sitter point as a power function of (1+z). Otherwise if (a — l) 2 — 4b < 0, the dark energy shows 
an oscillating behavior as 

, l=a , \ (a — l) 2 — 4b . c , , q , , 

y = y + 2B1 + z)~>~ cos ^ In 1 + z + 1 + z 3 . 36 

2 6a + b + o 



Now we write out the discriminant, 



Rip' V 

Combining the stability condition Eq. (|TT|) . 1 — t/j > and Rip' > we know, when 



(a I; 2 V, 1^- (-2 ^/]. (87) 



<7 2--5-f , (38) 



ip 4 V i?^' 

the dark energy has an oscillating behavior near the de-Sitter point. In this case, according to the definition in Eq. 
fllUl) . one has 



1 4B6., ,1=. 
w de = -l + - — (1 + z) 2 cos 
3yo 



V(q-l)^-4b i 

ln(l + z) + a 



y (3a + & + 6) 



(1 + z) 3 , (39) 



where a = arctan y/(a— l) 2 — 46/ (1 — a). The EOS parameter of dark energy also has an oscillating form. We have 
noticed that, on condition 1 — a > 0, the oscillating amplitude becomes smaller and smaller. Therefore, as time goes 
on, u>d e tends to be -1. If we take the e-folding number N — — ln(l + z) as the variable, the oscillating frequency is 



n*) = t/4--(2-^V- (40) 

IV. OSCILLATING IN THE MATTER ERA 

In the matter dominated era, z > 3 and R = 3if 2 , in order to match the local gravity tests, the following conditions 

/i-»i?-2A,Fi-M, (41) 
/a-M.-Fa-X^-H^'-X), 

must be satisfied. And the minimum point of the effective potential in the Eq. (IT91 is 

Rip - 2/i - f 2 T = 0. (42) 
So, we neglect the dark energy contribution, and we write out the expression of Ricci scalar 

R~ p = 3fh 2 (l + z) 3 . (43) 
In this case, according to Eq. (l23l) . the EOM in the matter era is 

d 2 y(z) , - 3 -^Wdy(z) , (2-V) 



dz 2 1 + z dz RiP'(l + z) 2 



y{z) (44) 



Using the method in the literature [2l[ , with a little modification, we solve the EOM near the minimum of the effective 
potential z = z + (z — z ), and z — z <C z . To first order in (z — z ), the Eq. (I44[) changes into 

, 2 + — 1 — + -— r^y z )=c z-z +rf l + z , (45) 

az^ 1 + zq dz (1 + zo)^ 
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where 



l - 4>(R ) 

a = - 3 -2R^ry (46) 

i?oV>'(i?o) ' [ ' 



3 , -4-^(i?o) + 5/i(i? )/i?o 3^"(i?o)(2-^(i? ) - fi(Ro)/Ro) 



(48) 



2 2R i)'{R ) 2^' 2 (i?o) 

6(F 2 (i? ) + ZRof^RoMiRp) - 6F 2 (R Q )3R (yj" (R ) - 2ft,'(R )) 

*P' 2 (R ) 

1 - jj(Ro) 6F 2 (Ro) WiCgoVgo 

2RotP>(Ro) ^'(iZo) 2R iP'(Ro) ' 1 j 



The solution of Eq. (J45J) is 



»(*) = f (1 + ^o) 3 + ^±^(, - , ) - |(1 + ,o) 3 + Cbe^^^, (50) 

where Co is a constant depending on initial conditions. Generally speaking, for different forms of tp , 1/R n for 
example, corresponds to (1 — ip)/Rip' ~ O(l) and e~ R corresponds to (1 — ip)/Rtp' ~ 0. Note that Rip' <C 1. 
Therefore, a ~ C(l), & ~ l/(Rip') ^> 1, c ~ d ~ (A/R)/(Rip r ) and ac/6 2 <C 1 which mean we can neglect the third 
term in Eq. (|50|) . In this case, the discriminant in the square root of Eq. (|50| must be negative, and the solution (|50|) 
becomes an oscillating form 

y(z) ~y Q + f(l + z ) 2 (2 ~ *>) + 2C e^7 CO s ($gg(* - z )) , (51) 

where yo — d(l + zq) 3 /b. Using Eq. (|13[) again, we give EOS parameter of dark energy near the redshift z — zq, 



c , „ \J\b - a 2 =» / \J\b - a 2 , ,\ , . 

u> de = -1 + ^—1 + z a 3 - Co^— sin ; (z - z ))■ 52 

3oy 3y y2(l + z ) y 

Since we care about the oscillating, we present the oscillating frequency, 

F(*o) * n T \ nr^ i (53) 

(1 + z )^/R ip' 

which corresponds to that in normal f(R) theory [2l|, where we have used the condition a 2 <C 46. 

If we hope the non-minimal coupling f(R) model improve the singularity problem, according to the expression (1531) .it 
is required that Rip'(l + z) must not be the decreasing function of Ricci scalar. In the matter dominated era, we 
ignore the normal fi(R) part. On the other hand, by considering think the scalar field ip stays at the minimal point 
of the effective potential, or equivalently R ~ p, we get 

Rtf ~ ~f%R 2 . (54) 

Recalling the stability condition (ITTl) . we get the constraint of f 2 

< -f^R 2 < 1. (55) 

On the other hand, to evade the singularity problem, considering 1 + z ~ R 1 ^ 3 , —f^R 7 ^ 3 must not be a decreasing 
function of Ricci scalar. Therefore we get the condition that improve the singularity, 

-iN- «> 



V. SPECIFIC COUPLINGS 



In the normal f(R) gravity, there is a singularity problem that the oscillating frequency is too large due to Rf"(R) <C 
1. We expect to reduce the oscillating frequency by adding a coupling term. Here, we consider some specific models. 
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A. logarithmic coupling 

Basing on the f(R) theory, we consider a logarithmic coupling 

f 2 (R) = l + a\n(R/c), (57) 

where c is the cosmological constant sharing the same unit with the cosmological constant A, and the parameter a 
is a small non-dimensional constant. In the matter era, for most of f(R) theory, to evade the solar system tests, the 
Rfi is set to be very small. Here, we set Rf" <C a <C 1 reasonably. Such a coupling has a virtue that it has little 
influence on the evolution of the early universe, for example, inflation and reheating era, as long as we take a proper 
coupling constant a. In the inflation and reheating era, considering c is a cosmological constant A, R/c ~ 10 121 , we 
get /a ~ 1 + af 21 In 10, especially when a < 1CP 4 , we get f 2 ~ 1 approximately. This recovers the normal f(R) theory. 
As for the local gravity, R/c ~ 10 5 , however, the detailed discussion, which we put in the Sec. VII, may be more 
complex. Here, we just set a to be a small constant without dimension. Therefore, the scalar ip is 

yj = f[- 2a| (58) 

where we have used the approximation R ~ p. According to the definition, ip' can be rewritten as 

So, we have the approximation ip ~ 1 — 2a and i?^>' ~ 2a. In this case, we can get the oscillating frequency in the 
matter era, 

1 lir 

F(zo) = —— — —i= = T , (60) 
(I + z )v2a J- 

which indicates that, the period T is proportional to the red-shift z if we consider a small logarithmic coupling. This 
solves the singularity problem in the matter era. 



B. power-law coupling 

We also consider the power-law coupling, assuming that 



MR) = !+(#] , Kill 



Rr, 



where R n is a constant with the same unit as M^. 



Firstly, we discuss the stability condition in power-law coupling, and we get 



= f"R - IPR^ 1 (£) ^ fi'R - 2n(n - 1) ( £) 



(62) 



In order to improve the large oscillating frequency, assuming (R/R n ) n 3> f['R in the matter era, according to Eq. 
([17]) . one has 

< -2n(n - 1) ( -^j < 1. (63) 

Then we have two condition 

0<n< 1, ( A ] < — i -. (64) 

2n(l-n) V ; 

Therefor, we do not consider the inverse power-law coupling n < which is used to mimic the dark matter in [42j], 
either the case n > 2. Using condition ([56l) . we get n < 13/3, which is always satisfied, if we consider stability 
condition. 

Now we set n = 1 — e in order that f 2 ^ 0, where e <C 1 is a small positive constant. Recalling the constraint 
arising from the non-minimal coupling scenario for reheating (43l \4l\ , 

A < lO-^ 4 , (65) 
Hi 

we know, in the matter era, though condition (jTTJ) and (|5T>|) is satisfied, (R/Ri) 1 ~ e is so small that it has no contribution 
to Rij) 1 . Hence, the power-law coupling is unavailable to improve the singularity problem. 
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C. exponential coupling 

Ones assume a power-law exponential coupling of the form 

/a=e(&)", (66) 
where the ( ) must be very small in the matter era to recover the normal f(R) theory. We get the Rip', 

R 



(n - 1) + n 

An 



(67) 



(n \ n— 2 r 
— 
it,, 

~tfi2-2n(n-l) 

Therefor the discussion is the same as the power-law coupling. 
There are also some other exponential coupling, such as 

h = 1 - e- K/c , (68) 

where the c is a constant which sharing the same unit and order with the cosmological constant A. In this model, 
fH ~ e~ R l c decreases rapidly that means Rip' is still very small in the matter era. So, the exponential coupling is 
unavailable to improve the singularity problem either. Using condition (|56l) . we get 



R/c, (69) 



which is bigger than 7/3. 



VI. NUMERICAL RESULT AND THE LATE TIME EVOLUTION 

In this section, using numerical calculation, we examine and certify the theoretical calculation about the oscillating 
frequency in the matter era. Then, using several possible coupling forms, we discuss the behavior EOS parameter in 
the de-Sitter universe. 

A. non-singularity with a small logarithmic coupling 

We set the coupling constant a = 10 -10 , and the initial condition z = 8, y(z ) = 2A, y'(z ) = . Therefore, in the 
matter era, seeing Eq. (|58[) and Eq. (|59[) , when the contribution of /i is insignificant comparing with the coupling 
term, the coupling constant will dominate. This will suppress the oscillating in the early universe. Fig. 1 shows 
the EOS parameter of dark energy as a function of red-shift. In the late time, with a small coupling, we can not 
distinguish the result from that in normal f(R) theory without coupling. However, in the matter era, there are some 
differences. Fig. 2 shows the oscillating at different red-shift. Note that, as the red-shift gets bigger, the oscillating 
period getting bigger. It reads as: 
redshift numerical calculating period theoretic predict 

z=4 T= 0.00042 0.00044 

z=5 T= 0.00053 0.00053 

z=6 T= 0.00062 0.00062 

z=7 T= 0.00071 0.00071 



Comparing with the result in the literature 2l|, there is no more singularity in this model. Form the Table. 1, we 
get the relationship 

AT - Az, (70) 



which can be also got from the theoretical result Eq. (|60|). Note that, when 2 = 4, there is some differences between 
the numerical and the theoretical calculation. This is because our approximation ip — > 1 and Rip' < 1 is no longer 
valid in the low redshift region. However when z > 4, our theoretical calculation is still valid. 



h=2, a=Kr 10 , z =8, t.=R-\ c(1-e~ R,c ), f =1+a ln(R/c) 

-0.94 1 1 1 1 1 1 1 r 
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FIG. 2: The EOS parameter of dark energy at different red-shift, 
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N 

FIG. 3: The EOS parameter of dark energy in exponential model with a exponential coupling, where the parameter A = 2 

B. late time evolution in toy models 

If we ignore the singularity problem in the early universe and care more about the late time universe, we want to 
check how the non-minimal coupling wound influence the late time universe. To match the solar system tests, when 
Ricci scalar is large, f2(R) must tend to be 1 while considering the coupling form. We still consider a exponential 
model, and with a exponential coupling: 

A(R) = R-\c(l-e- R / c ), (71) 
h(R) = l-e- R / c . (72) 

When A > 1, in the matter era, e~ R l c <C 1 which assure /a — > 1. Fig. 3 shows the detail, when — 1 < N < 3, we can 
distinguish the results from that in normal f(R) theory. When N > 3, the result is nearly the same. In this model, 
we set A = 2 and the initial red-shift is zo = 3. 
We also consider the model 

f 1 (R) = (R-\c)e x ^" (73) 

proposed in [l8| recently. Therefore we add new exponential coupling form 

f 2 (R) = e Al ^)" (74) 

when A > 1, in the matter era, (c/R) n — > that assure f% — > 1. where the parameters we set is A = 1, n = 4 and 
Ai = 0.5. the initial red-shift is zq = 12. Fig. 4 shows the details. In this coupling model, the oscillating amplitude 
is reduced comparing with the normal f(R) theory. 

VII. LOCAL TESTS 

In this section, we consider local tests of the non-minimal coupling f(R) model. We will use the similar method by 
Hu and Sawicki [14j. We consider a spherically symmetric isotropic metric, 

ds 2 = - [1 - 2A(r) + 2B{r)\ dt 2 + [1 + 2A(r)} (dr 2 + r 2 dil), (75) 

where we assume A -C 1 and B <C 1 near the source, especially, in the GR limit, B — > since the solar system tests 
provide the strongest limits on B. 
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FIG. 4: The EOS parameter of dark energy in new exponential model with a new exponential coupling 



Generally, according to the definition of the Ricci tensor, we get 

R° = W 2 (A-B), (76) 
R=-2V 2 (A + B). (77) 

And form the modified Einstein equation ([2]), we get the time-time component of the field equation 

i>Rl-\h{R) + U^ + dfy = -f 2 K 2 p. (78) 

If we consider the static solution, combining this equation with the trace equation (|14j) . we get the time-time component 
of Ricci tensor 

3R° = J fl + * hK2p + R. (79) 



Then we get 

fi + 4/ 2 ^ 2 p 



V 2 A 


- V 






1 

12 


(* 


V 2 B = - 




^ = 


1 

12 


(- 


5R 



(80) 
(81) 



In the high curvature region, for general f(R) model, we have f±—>R and f[ — > 1. On the other hand, taking proper 
form of /2, which has the limit /2 — > 1 and f 2 in the high curvature region, we have ij) — > 1. So we consider 

VM«-i-^ 2 p, (82) 

V 2 £«i( K 2 P -i?). (83) 



In the same limit, from the trace equation of the field equation (|14[) . we have 



V 2 V« 1(R-k 2 p ). (84) 



Therefore, ^ together with B has a solution 



i>(R)+B{r)= ai + ^. (85) 
r 
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We assume the solution remains finite at r — 0, a 2 = 0, and when r — > 00 we assume B(oo) 0. So, we have 

B(r Q ) = -[ip(R Q )-iP(Roo)] (86) 

It is difficult to find an accurate solution of the Ricci scalar. However, we can consider the field stays at the minimal 
point = n 2 p approximately, and the deviation from the minimum is 

i?« = 3D^(i2) + {i>- l)R - 2(/i -R) + (/ 2 - 1)k 2 P , (87) 

where the last three terms we know is small, so we care the field gradients V 2 V>- And we define R = RW+RW. A 
sufficient condition for R^ <C R^ is the compton condition 

ip'\ R = K * p dfp < p. (88) 

If this condition is always satisfied, we consider R « n 2 p. If this condition is violated beyond some outer radius, such 
as the sun radius, this approximation is not valid until the compton condition is locally satisfied again outside the 
sun. Then the the field stays at the minimum again. 

Anyway, the curvature near the sun is high, so we think ip(R Q ) w 1 comparing with ip(po). On the other hand, we 
use the galaxy radius instead of 00, considering the field stay at the minimum point Roo = R g = n 2 p g - So we get 

B(r ) « i>{R a ) - 1 (89) 
Finally the deviation from the GR metric is given by 

7-1^ (90) 

where A we take is the Newtonian potential of the sun $jv = GM/r = 2.12 x 1CP 6 . The tightest experimental bound 
on the PNP parameter is given by I7 — 1| < 2.3 x 1CP 5 [48l450l |. Therefore we get the constraint 

rl>(Rg) - 1 = (f[(R g ) - 1) - 2f 2 p < 0.5 x lO" 10 . (91) 

Generally, (f[(R g ) — 1) < 0.5 x 10~ 10 is the constraint of normal f(R) component. If we require 

V2{R a )p a < 0-5 x 10- 10 , (92) 

the coupling model may evade the solar system test. For logarithmic coupling, f 2 = 2ap g /R g w 2a, the constraint is 
a < 1Q- 1Q . 



VIII. CONCLUSIONS AND DISCUSSION 



We consider a non-minimal coupling between geometry and matter in the f(R) theory. We find this theory also has 
a independent scalar ip = f[ + 2Lf 2 similar to the scalar /' in the normal f(R) theory. We give its stability condition 
(fT7|) . which is also similar to normal f(R) theory. 

We discuss the evolution of dark energy in the de-Sitter universe, when the condition (p?81 is satisfied, the dark 
energy density will have oscillating behavior. Meanwhile, we find the dark energy is always oscillating in the matter 
era, with an oscillating frequency F(z) cx 1/(1 + z)\fRip' . Adding the coupling part into the f(R) theory, under the 
condition that — ^„ fl < |, the oscillation is suppressed successfully by the coupling term in the matter era. 

Three specific coupling forms have been studied, and it turns out that the logarithmic coupling, with a oscillating 
period AT ~ Az, is a good choice. The power law coupling also can suppress the oscillating becoming strong. 
However, the reheating constraint < 10 -104 gives a overlarge energy scale. And this will disappointedly lead to 
very small R/R\, which makes it difficult to reduce the oscillating frequency in the matter era. The exponential 
coupling /a = 1 — e~ R can not satisfy the increasing function condition (|69l) . 

Through numerical calculation, we verified our theoretical calculation of the oscillating frequency in the matter 
era, and found only when z < 4 the numerical calculation has some differences with the theoretical calculation. This 
is because our approximation ip — > 1 and Rip' < 1 is no longer valid in the low redshift region. However when 
z > 4, our theoretical calculation is still valid. We take two toy models as examples, and show that the non-minimal 
coupling does not change the cosmic evolution of f(R) model in the late universe. Especially when Ina > 4, we can 
not distinguish between coupling model and the f(R) model. 

In the end, we discuss the local test constraint on the non-minimal f(R) model. Generally, we obtain the constraint 
ip g — l< 10~ 10 . in the logarithmic coupling model, this constraint forces us to set the coupling parameter a to 10~ 10 , 
which leads a small oscillating period T ~ 10 -4 , seeing Fig. 1. However, its oscillating frequency is decreasing as the 
redshift increasing and the amplitude is not too large comparing with the normal f(R) theories. So, the singularity 
problem in normal f(R) model is no longer so serious when we consider the non-minimal coupling. 
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